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Direct computer simulations of electromagnetic scattering by discrete random media have 
become an active area of research. In this progress review, we summarize and analyze our 
main results obtained by means of numerically exact computer solutions of the macro- 
scopic Maxwell equations. We consider finite scattering volumes with size parameters in 
the range [20,60], composed of varying numbers of randomly distributed particles with 
different refractive indices. The main objective of our analysis is to examine whether all 
backscattering effects predicted by the low-density theory of coherent backscattering (CB) 
also take place in the case of densely packed media. Based on our extensive numerical data 
we arrive at the following conclusions: (i) all backscattering effects predicted by the 
asymptotic theory of CB can also take place in the case of densely packed media; (ii) in the 
case of very large particle packing density, scattering characteristics of discrete random 
media can exhibit behavior not predicted by the low-density theories of CB and radiative 
transfer; (iii) increasing the absorptivity of the constituent particles can either enhance or 
suppress typical manifestations of CB depending on the particle packing density and the 
real part of the refractive index. Our numerical data strongly suggest that spectacular 
backscattering effects identified in laboratory experiments and observed for a class of high- 
albedo Solar System objects are caused by CB. 

Published by Elsevier Ltd. 
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1. Introduction 

Two optical phenomena, viz., the sharp and narrow 
brightness opposition effect (BOE) and the polarization 
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opposition effect (POE, otherwise known as the azimuthal 
asymmetry of the polarized backscattering cone), have 
been observed for a class of high-albedo Solar System 
bodies (see, e.g., Refs. [1-4] and references therein) as 
well as identified in laboratory measurements for parti- 
culate media with large values of the volume packing 
density [5-10]. The BOE manifests itself as a spike-like 
intensity peak centered at exactly the backscattering 
direction, while the POE is observed in the form of a 
sharp asymmetric negative-polarization feature with a 
minimum at a phase angle comparable to the angular 
semi-width of the BOE. Both features have been attrib- 
uted to the effect of coherent backscattering (CB) in 
particulate media (see, e.g., Refs. [4,11] and references 
therein). It should be kept in mind, however, that the 
interference concept of CB is strictly applicable only in the 
asymptotic limit of infinitesimally small packing density 
[11-14]. Therefore, one faces a fundamental question of 
whether the effect of CB can also account for the specific 
backscattering phenomena observed in the case of den- 
sely packed media. 

A definitive answer to this question can only be 
obtained using extensive numerically exact computations 
of light scattering characteristics of media composed of 
varying numbers of randomly positioned particles [13,14]. 
The steadily increasing power of computers and the 
availability of efficient numerical techniques have 
recently facilitated the emergence of an accurate quanti- 
tative approach to solving this complex problem based on 
direct computer solutions of the macroscopic Maxwell 
equations [15-43], This approach can be used, in parti- 
cular, to assess predictions and conditions of applicability 
of different approximate theories and idealized physical 
concepts [13,14,32,38]. 

In this progress review, we analyze the behavior of 
backscattering characteristics of macroscopic volumes of 
discrete random medium with varying packing density 
using numerically exact data obtained recently by directly 
solving the macroscopic Maxwell equations. The main 
objective is to examine whether the optical effects pre- 
dicted by the low-density theory of CB also take place in 
the case of densely packed media. We provide a summary 
of the main results obtained in our recent publications as 
well as discuss new inferences based on numerical results 
for particulate volumes with very high packing densities 
or composed of absorbing constituent particles. 

2. Theory 

The effect of CB is a fascinating optical phenomenon, 
which can survive essentially any degree of randomness 
of particle positions [11,12]. The interference origin of CB 
has its roots in the far-field version of the Foldy-Lax 
equations [11] and is explained schematically in Fig. la, 
which shows a random multi-particle group illuminated 
by a plane electromagnetic wave incident in the direction 
of the unit vector n iB and observed from a very large 
distance. If the direction of observation, given by the unit 
vector n ote , is far from the exact backscattering direction 
defined by the unit vector -n iu then the average effect of 
the interference of a pair of conjugate scattered waves 


going through a chain of n particles in opposite directions 
is zero, owing to the randomness of particle positions. 
Consequently, the detector of electromagnetic energy 
measures some average, incoherent (or diffuse) intensity. 
However, when the phase angle a (i.e., the angle between 
the vectors n oi)s and -n iu ) is zero, the phase difference 
between the conjugate paths involving any chain of 
particles vanishes, and the interference is always 
constructive. 

The most well-known manifestation of the interfer- 
ence origin of CB is the BOE in the form of a narrow peak 
of intensity centered at oe=0 r (see, e.g., Refs. [6-8]). On the 
other hand, the physical explanation of the interference 
origin of the POE is more complicated [44-46] and is 
given in Fig. lb. 1 Particles 1-4 lie in a plane normal to the 
illumination direction and are assumed to have sizes 
much smaller than the wavelength. Particles 1 and 2 lie 
in the scattering plane (i.e., the plane containing both unit 
vectors n iB and n ote , shown by yellow), while the line 
through particles 3 and 4 is normal to this plane. If the 
incident light is unpolarized then both magenta trajec- 
tories cause scattered light polarized negatively with 
respect to the scattering plane, whereas both blue trajec- 
tories cause positively polarized scattered light. The phase 
difference between the conjugate magenta trajectories is 
always zero, while that between the blue trajectories is 
zero if a=0° and oscillates rapidly with increasing a. 
Therefore, on average, CB enhances the contribution of 
the negatively polarizing scattering trajectories over a 
wider range of phase angles than that of the positively 
polarizing trajectories. The result is a negative-polariza- 
tion minimum at a small phase angle comparable to the 
angular semi-width of the coherent BOE [44,49]. The fact 
that only specific particle configurations contribute to the 
POE often makes the latter less pronounced than the BOE. 

Fig. 2 shows the results of computations performed in 
[49] by using the asymptotic low-density theory of CB in 
the case of a semi-infinite homogeneous slab composed of 
sparsely distributed non-absorbing Rayleigh scatterers 
[50,51]. One can clearly see the most obvious manifesta- 
tions of the effect of CB such as a narrow peak of intensity 
centered at zero phase angle and a sharp asymmetric 
minimum of negative polarization at a small phase angle 
comparable to the angular semi-width of the coherent 
intensity peak. 

The interference concept of CB explicitly relies on 
assigning a phase to the wave scattered by a particle. 
This implies that particles must be located at large 
distances from each other, in the so-called far-field zone 
(see, e.g., Refs. [11,12]). However, as we have stressed in 
Section 1, many laboratory measurements for densely 
packed particle suspensions and particulate surfaces as 


1 It is worth clarifying that this optical phenomenon was introduced 
in Refs. [47,48] as a way of explaining the wide, nearly parabolic 
negative branch in the polarization phase curves exhibited by the 
majority of atmosphereless Solar System bodies [4], However, it was 
shown in Refs. [44,49] that CB results in the highly asymmetric spike- 
like POE rather than in the much more ubiquitous, nearly symmetric 
negative-polarization branch (NPB). The very term "POE” was intro- 
duced in Ref. [44], 
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Fig. 1 . (a) Schematic demonstration of the interference origin of the CB effect, (b) Schematic explanation of the CB origin of the POE. (c) Scattering by a 
macroscopic spherical volume randomly filled with small spherical particles, (d) Polarization of scattered light for a single ice sphere with different 
values of the size parameter kr. (e) Polarization of scattered light for a single sphere with several combinations of the refractive index m and size 
parameter kr. 


well as observations performed for a class of high-albedo 
Solar System bodies appear to exhibit the unique features 
predicted by the low-density theory of CB. Therefore, it is 
necessary to perform numerically exact computations of 
light scattering based on a direct computer solution of the 
macroscopic Maxwell equations in order to examine the 
applicability of the sparse-medium concept of CB to 
densely packed media in which particles are often in 


direct contact with each other instead of being separated 
by large distances. 

3. Numerically exact computer modeling 

Our model of random particulate medium is a sphe- 
rical volume of radius R filled with N identical non- 
overlapping spherical particles of radius r (see Fig. lc). 
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Fig. 2. Angular profiles of the enhancement factor (black curve) and the degree of linear polarization (gray curve) for a semi-infinite homogeneous slab 
composed of sparsely distributed non-absorbing Rayleigh scatters and illuminated by normally incident unpolarized light [49], The enhancement factor 
is defined as the ratio of the total scattered intensity to that of the incoherent (diffuse) background. The angular parameter q is defined as the product of 
the phase angle, the wave number, and the mean free path of light in the scattering medium. 


The size parameter of the volume is kR, and the particle 
size parameter is kr, where k is the wave number. To 
model statistically random and uniform particle positions 
within the spherical volume consistent with the funda- 
mental assumption of ergodicity [11], we follow the 
approach pioneered in Ref. [22]. Specifically, we use one 
realization of an N-particle group generated randomly 
according to the procedure described in [52], and then 
average all optical observables of interest over the uni- 
form orientation distribution of this configuration with 
respect to the laboratory coordinate system. Although the 
mutual positions of the N particles with respect to each 
other remain the same, they are sufficiently “random” to 
begin with. Therefore, this simple approach yields, in 
effect, an infinite continuous set of random realizations 
of the scattering volume while enabling one to use 
the efficient orientation averaging technique afforded by 
the superposition T-matrix method (STMM) [53-55]. The 
STMM represents a direct computer solver of the macro- 
scopic Maxwell equations for an arbitrary multisphere 
configuration [56]. Within the range of its numerical 
convergence, the corresponding public-domain T-matrix 
computer code [57] yields results with a guaranteed 
accuracy, which makes it numerically exact. 

It is obvious that by its very construct, this scattering 
model cannot be expected to reproduce exactly the infinite 
diversity of morphologies of particulate media encountered 
in laboratory and natural conditions. However, it proves to 
be sufficiently representative to permit a robust and instruc- 
tive analysis of the effects of packing density on multiple 
scattering of light by discrete random media. 

It is further assumed that the statistically random 
particulate volume is illuminated by a parallel quasi- 
monochromatic beam of light propagating in the direction 
of the unit vector n i(/ (Fig. lc). The observer is located in 
the far-field zone of the entire spherical volume in the 
direction of the unit vector n obs . Since all scattering and 
absorption properties of the particulate volume are aver- 
aged uniformly over all orientations of the N-particle 
group, they depend only on the phase angle a provided 
that the scattering plane is used for defining the Stokes 
parameters of the incident and scattered light. The far- 


field transformation of the Stokes parameters /, Q, U, and V 
upon scattering by the entire particulate volume is then 
formulated in terms of the 4x4 real-valued Stokes 
scattering matrix [11,56,58] 


jsca - 


a i(a) 

b i(a) 

0 

0 " 


- jinc - 

Qsca 

oc 

bi(a) 

a 2 (a) 

0 

0 


Qinc 

[jsca 

0 

0 

a 3 (a) 

b 2 (a) 


U inc 

ysca 


0 

0 

-b 2 (a) 

a 4 ( a) 


v inc . 


In theory [56,59], the off-block diagonal elements of 
the scattering matrix would be zero only if the primordial 
N-sphere configuration had a plane of symmetry. We 
found, however, that in all the cases considered below 
the magnitude of the minor matrix elements denoted in 
Eq. (1) by zeros is much smaller than that of the major 
elements. In other words, the zeros denote the scattering 
matrix elements negligibly small (in the absolute sense) 
relative to the other elements at the same scattering 
angles. The (1, 1) element a^a), called the phase function, 
is normalized according to the standard integral condition 

1 f n 

-/ a^(a)smadc/. = 1. (2) 

2 Jo 


The elements of the scattering matrix can be used to 
define conventional optical observables corresponding to 
different types of polarization state of the incident radia- 
tion used in remote-sensing, in situ, and laboratory 
particle characterization (see, e.g., Refs. [60-66] and 
references therein). Specifically, if the incident light is 
unpolarized then the phase function characterizes the 
angular distribution of the far-field scattered intensity, 
while the ratio -tq/cq gives the corresponding degree of 
linear polarization. If the incident radiation is polarized 
linearly in the scattering plane (i.e., Q mc =I mc , U inc =V inc =0) 
then the linear polarization ratio /. i L is defined as the ratio 
of the cross-polarized and co-polarized components of the 
scattered intensity 


ai(«)— q 2 («) 

L " a^(a)+2b^(a)+a 2 (a) ' 


( 3 ) 
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If the incident radiation is polarized circularly in the 
counterclockwise direction when looking in the direction 
of propagation (i.e., V inc =I inc , Q‘ nc =U inc =0), then the circu- 
lar polarization ratio ji c is defined as the ratio of the 
same-helicity and opposite-helicity components of the 
scattered intensity 


fi c (a) = 


a 1 (a) + a4(a) 
ai(a)-a 4 (a) ' 


(4) 


By using the numerically exact computer solver of the 
macroscopic Maxwell equations, we have performed 
extensive computations of the scattering characteristics 
of particulate volumes with size parameters kR in the 
range from 20 up to 60, filled with non-absorbing 
spherical particles with real-valued refractive indices 
m = 1.194 (representing latex in water), 1.31 (represent- 
ing liquid water and water ice at visible wavelengths), 
and 1.4, 1.5, 1.6 (representing different mineral sub- 
stances) [33,34]. Depending on kR, the number of mono- 
mers N varies from 1 up to 900 (for m=1.31 and kR= 30). 
Also we present some numerical data for absorbing 
particles with m R =Re(m) = 1.31 and 1.5 and m,=Im(m)= 

0.01, 0.1, and 0.3. 

The corresponding values of the packing density 
depend on the way the particle volume fraction is 
defined. If R is the radius of the sphere enclosing all 
constituent particles completely then the packing den- 
sity p is defined as the ratio of the cumulative volume of 
the constituent particles to the total enclosed volume: 
p = Nr 3 /R 3 . However, this definition underestimates the 
actual packing density inside the scattering medium 
since the constituent particles are not allowed to cross 
the bounding sphere (Fig. lc). Therefore, in what follows 
we define the packing density as p = Nr 3 /(R-r) 3 , where 
R—r is the radius of the sphere enclosing all constituent- 
particle centers. 

It should be kept in mind that a direct computer 
modeling of the POE can be quite problematic because 
the degree of linear polarization of light singly scattered 
by a constituent particle (for unpolarized incident light) 
often has a negative branch at small phase angles [67,68]. 
As an example, Fig. Id shows the computed values of 
polarization for an ice sphere with the refractive index 
m = 1.31 and several values of the size parameter kr. We 
see that in all cases, the computed polarization has a 
negative branch in the range of phase angles [O', 20°]. This 
fact makes it non-trivial to distinguish between the singly 
and multiply scattered negative-polarization contribu- 
tions in a direct numerical solution of the Maxwell 
equations. In order to eliminate this uncertainty, for each 
value of the refractive index studied we have selected a 
size parameter value yielding a single-scattering polariza- 
tion contribution with a wide horizontal “shelf’ of near- 
zero values at small phase angles. The resulting {refractive 
index, size parameter} pairs and the corresponding single- 
scattering polarization curves are shown in Fig. le. One 
can clearly identify indeed a horizontal shelf of near-zero 
values of polarization in the range of phase angles 
0° < a < 30°, which makes any multiple-scattering polar- 
ization contribution easily discernable and quantifiable. 
Note that this approach was proposed in our recent 


publications [33,34] and was used to obtain the first ever 
numerically exact display of the entire suite of back- 
scattering optical effects implied by the asymptotic low- 
density theory of CB. 

4. Numerical results and discussion 

In this section, we present and analyze numerical 
results obtained for both non-absorbing and absorbing 
constituent particles. 

4.1. Scattering volumes filled with non-absorbing particles 

Figs. 3 and 4 provide a representative subset of our 
extensive numerical results computed for non-absorbing 
constituent particles. In some cases, we also show the 
results computed for a single isolated spherical particle. 
Since we are interested in the backscattering behavior of 
all optical observables, we consider only the range of 
phase angles 0° < a < 30°. 

Fig. 3a depicts the results of computations for m = 1.31, 
kR= 40, kr= 2, and packing densities p varying between 
1.5% (N=100) and 11.7% (N= 800). For reference, we also 
show the corresponding singe-particle data. We see that 
these numerically exact results are in perfect qualitative 
agreement with the predictions of the low-density theory 
of CB. Indeed, 

1. The normalized scattered intensity afia)lafi0) exhibits 
backscattering peaks rapidly developing with N. The 
angular widths of these peaks are approximately 
the same. 

2. In the phase-angle range 0° < a < 30°, the degree of 
linear polarization ~bfia)lafiot) is equal to zero for 
N=1 but rapidly develops a pronounced minimum 
with growing N caused by the increasing amount of 
multiple scattering. The phase angle of minimal polar- 
ization, a min , is virtually independent on N and is 
comparable to the angular semi-width of the coherent 
BOE. Furthermore, the angular shape of the negative- 
polarization minima is asymmetric, with a min being 
significantly closer to zero than to the phase angle at 
which polarization switches sign from negative to 
positive (the so-called inversion angle). As such, the 
polarization minima in Fig. 3a have the same basic 
morphology as the POE predicted by the low-density 
theory of CB, although the polarization curve in Fig. 2 
is based on the assumption of a phase-angle-indepen- 
dent ladder contribution and, as a consequence, shows 
no inversion. 

3. The backscattering peaks in the linear and circular 
polarization-ratio curves are absent in the case of a 
single particle but develop and rapidly grow with 
increasing N. The angular widths of all these peaks 
are approximately the same and independent of N. 

In order to analyze the dependence of these optical 
observables on the volume size parameter kR, we present 
in Fig. 3b the computational data obtained for m=1.31, 
N= 800, and kR= 30, 40, and 60. It is seen that the angular 
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N = 100 
N = 200 
N = 300 
N = 400 


N= 500 
N= 600 
N= 800 
N = 1 


b 



kR = 30, N = 800 
kR = 40, A/ = 800 


kR = 60, N = 800 
kR = 30, N= 100 


C 



Phase angle (deg) Phase angle (deg) Phase angle (deg) Phase angle (deg) 


m R = 1.31 
m R = 1.4 


m R =1.5 
m R =1.6 


Fig. 3. (a) Backscattering characteristics of a kR = 40 scattering volume randomly filled with a number N of kr= 2 spherical particles. The packing density 
varies between p = 1.5% (N= 100) and p = 11.7% (N=800). The particle refractive index is fixed at m = 1.31. (b) Backscattering characteristics of 
particulate volumes with size parameters in the range 30 <kR< 60. The particle refractive index is fixed at m= 1.31 and their size parameter is fixed at 
kr= 2. (c) Backscattering characteristics of a particulate volume with kR=40 and N=500. The particle refractive index m varies from 1.31 to 1.6. The 
respective values of the particle size parameter are listed in the legend of Fig. le. The packing density p varies between 7.3% (m= 1.31) and 4.1% (m = 1.6). 


widths of all the backscattering peaks and the angles of 
minimal polarization decrease with increasing kR 
approximately as 1/fcR. 


Also shown for comparison are the results for kR = 30 
and N= 100. In this latter case, the particle packing density 
p is approximately the same as in the case of kR = 60 and 
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Fig. 4. Optical characteristics of kR=20 (panels a and b) and 40 (panel c) scattering volumes randomly filled with different numbers N of spherical 
particles with (a) m=1.31, (b) m-1.5, and (c) m = 1.194. The packing density p in panel (c) varies between 5.8% ( /V = 5 0 ) and 46.7% ( /V= 400). 


N= 800 (3.6% and 3.3%, respectively). It is seen that for the 
same volume fraction, the backscattering enhancement in 
the intensity and polarization ratios as well as the depth 
of the negative-polarization branch increase with increas- 
ing kR, while their angular widths decrease approximately 
by a factor of two. 


In Fig. 3c, we show the results for the volume size 
parameter kR= 40, the number of monomers N= 500, and 
refractive indices m = 1.31, 1.4, 1.5, and 1.6 (the corre- 
sponding values of the particle size parameter kr are given 
in Fig. le). One can see very similar opposition effects for 
all four refractive-index values, which illustrates the 
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universal interference nature of CB. Quite remarkably, the 
angular widths of all the backscattering peaks are vir- 
tually independent of the particle refractive index, while 
the negative-polarization features are almost identical. 

From Fig. 3, we can conclude that all curves of normal- 
ized scattered intensity a 1 (a)/a 1 (0) and linear polarization 
— f>i(a)/ai(a) for N> 1 reveal the BOE and the POE, respec- 
tively. Indeed, all traits in the behavior of the normalized 
intensity peaks and the polarization minima testify to their 
interference multiple-scattering origin. Specifically, the 
intensity peaks and polarization minima are absent in the 
single-particle curves but are present when N> 1. The 
angular semi-width of the normalized intensity peak is 
comparable to the phase angle of minimal polarization 
a min . Both of them are independent of N and m and are 
inversely proportional to the value of the volume size 
parameter kR. Furthermore, the polarization minima 
become increasingly asymmetric with increasing kR and 
approach the classical angular profile shown by the gray 
curve in Fig. 2. Note, that the linear and circular polariza- 
tion-ratio curves are even more indicative of their multiple- 
scattering origin since they are absent completely for a 
single spherical particle. For a fixed kR, the angular widths 
of the backscattering peaks in the p L and fi c curves are 
approximately the same and are independent of the number 
of particles and their refractive index. Furthermore, they 
scale as 1/fcR, thereby corroborating their CB nature. 

The effect of increasing the number of particles N in a 
volume can be expected to be twofold. Initially it stimu- 
lates multiple scattering and thus enhances the classical 
manifestations of diffuse radiative transfer (RT) and CB. 
Eventually, however, it can cause changes in the scatter- 
ing patterns not implied by the low-packing-density 
theories of RT and CB. One should, therefore, expect that 
the RT and CB interpretation of the numerically exact 
T-matrix results must become inadequate when p becomes 
sufficiently large. 

In order to analyze the effects of increasing N, we 
present the STMM results for kR= 20, m = 1.31, kr= 2, and N 
varying from 50 (p ~ 7%) up to 345 (p ~ 47%) (Fig. 4a), as 
well as for m R = 1.5, kr= 1.76, and N varying from 100 
(p ~ 9%) up to 530 (p~48%) (Fig. 4b). We see that, 
indeed, starting from N=300 (p ~ 41%) the curves in 
Fig. 4a develop high-frequency interference ripples typi- 
cal of a single spherical particle with a size substantially 
greater than the wavelength [56]. Similar behavior is seen 
in Fig. 4b starting from N=400 (p ~35%). Obviously, this 
behavior is not predicted by the low-density theories of 
RT and CB. Nevertheless, the results of our computations 
do demonstrate that the predictions of these asymptotic 
theories can survive substantial volume packing densities 
typical of particulate surfaces and particle suspensions. 

It should be noted that the dense-packing behavior 
of the backscattering characteristics considered is not 
universal, as the data shown in Fig. 4c demonstrate. 
Indeed, the STMM results for kR= 40, m = 1.194, kr= 3.81 
(cf. Fig. le), and N varying from 50 (p ~ 6%) up to 400 
(p ~ 47%) exhibit the “standard” CB behavior despite 
packing densities reaching 47%. Apparently, this behavior 
is somehow related to the optical “softness" of the 
constituent particles. 


4.2. Scattering media composed of absorbing particles 

It is well known that absorption can diminish various 
manifestations of multiple scattering significantly (see, 
e.g., [11,25]). Therefore, it is important to verify whether 
the CB multiple-scattering interpretation of the back- 
scattering features identified for scattering media com- 
posed of non-absorbing particles is consistent with the 
way in which these features should be expected to change 
upon increasing absorption. Note that earlier in [25], 
using the results of computations for m R = 1.32 and 
p = 22%, it was concluded that increasing absorption 
suppresses such optical effects as BOE and depolarization. 

Some of the results of our extensive computations are 
presented in Fig. 5. The computations have been performed 
for a kR = 30 scattering volume composed of spherical 
particles with m R =1.31 (Fig. 5a), 1.5 (Fig. 5b), and m,= 0, 

0. 01. 0.1, and 0.3. The values of the particle size parameter 
were adopted to be the same as in the case of m/= 0 (i.e., 
kr= 2 for m R =1.31 and 1.76 for m R = 1.5). In the case of 
m R = 1 .3 1 , the number of constituent particles N was varied 
from 1 00 up to 400, which corresponds to the variation of 
the packing density p from 3.6% up to 14.6%; for m R = 1.5, 
N was varied from 100 (p = 2.4%) up to 600 (p = 14.5%). 

Fig. 5a presents the results for the real part of the 
refractive index m R =1.31. We see that: 

1. In the range of phase angles 0°<a<10 D , the ampli- 
tude of the backscattering peak in the normalized 
scattered intensity a 1 (a)/ai(0) caused by CB decreases 
with increasing absorption. 

2. The depth of the negative branch in the curves of linear 
polarization -h 1 (a)/a 1 (a) decreases monotonously with 
increasing m , ; when m/= 0.3, the polarization curve 
changes its shape. 

3. Increasing absorption diminishes progressively the 
values of the linear and circular polarization ratios; 
when m,= 0.3, both ratios become very small. 

All these results can easily be interpreted in terms of the 
decreasing contribution of multiple scattering with increas- 
ing absorption. It is known that absorption terminates long 
scattering paths and thereby can be expected to reduce the 
range of phase angles affected by CB [69]. This can be seen 
indeed in our results which exhibit, in particular, a decrease 
in the width of the BOE and the phase angle of minimal 
polarization a min with decreasing m h 

Somewhat different results are obtained in the case of 
the real part of the refractive index m R = 1.5 and the number 
of constituent particles N=100 (Fig. 5b). In this case: 

1. Increasing the imaginary part m, from 0 up to 0.1 
results in a stronger backscattering intensity peak 
a 1 (a)/a 1 (0). Furthermore, the backscattering enhance- 
ment is still quite pronounced when m;=0.3, which is 
indicative of a significant residual contribution of 
multiple scattering even in the case of such strong 
absorption. 

2. The behavior of the depth of the negative-polarization 
branch with increasing absorption is not monotonous: 
it first increases and then decreases. 
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a m R = 1.31, kR = 30 



b m R = 1.5, kR = 30 



Fig. 5. Scattering characteristics of a /<K=30 scattering volume randomly filled with N identical spherical particles, (a) m K =1.31, N = 1 00 (upper row), and 
N=400 (lower row); (b) m R =1.5, N= 1 00 (upper row), and N=600 (lower raw). The imaginary part of the refractive index varies in the range 0 < mi <0.3 
according to the color legend in the top right-hand diagram. 


3 The polarization ratios /.i L and fi c also demonstrate a 
significant growth as m , increases from 0 up to 0.1, 
followed by a precipitous decrease. At the same time 
the widths of the backscattering peaks hardly change. 

However, in the case of m R = 1.5 and N=600, we see a 
more monotonous decrease of all the manifestations of 
the effect of CB with increasing absorption. 

Thus, we have to conclude that increasing absorption 
can either enhance or suppress manifestations of the CB 
effect depending on the particle packing density and the 
real part of the refractive index. This is not surprising, 
since increasing m/ changes not only the absorption 


properties of the constituent particles but also their 
scattering (including polarization) properties. 

5. Concluding remarks 

It must be emphasized once again that our model of a 
multi-particle scattering volume cannot be expected to 
reproduce exactly the diverse complex morphologies of 
particulate media encountered in laboratory and natural 
conditions. As noted above, the interference base for a 
finite scattering volume is controlled by its size parameter 
kR, whereas that for an optically thick, non-absorbing or 
weakly absorbing layer is controlled by the transport 
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mean free path [12], The latter can be much greater than 
the kR values used in our computations and thereby yield 
much narrower opposition effects. However, our model 
has proved to be remarkably robust and permits a 
thorough qualitative and semi-quantitative analysis of 
the effects of packing density on multiple light scattering 
in discrete random media. The critical advantages of the 
direct computer solution of the Maxwell equations is that 
it yields numerically exact results, does not involve the 
simplifying yet questionable assumption of a small pack- 
ing density, and allows the researcher to vary all physical 
parameters of the scattering medium one a time [13,14]. 
This allows one to trace the onset of multiple scattering 
as the particle number N deviates from one, the evolution 
of the various manifestations of CB with increasing N, and 
the eventual onset of packing-density effects [70,71] 
distorting substantially the angular profiles of the BOE, 
POE, and polarization ratios with respect to the profiles 
predicted by the asymptotic low-density theory of CB. 

On the basis of the results of our extensive numerically 
exact computer simulations [20,22,25,31,33,34], we arrive 
at the following conclusions: 

1. All backscattering effects predicted by the low-density 
theory of CB can also take place in the case of a densely 
packed medium. The relatively small optical diameters 
of the particulate volumes considered cause the angu- 
lar widths of the BOE and POE to scale as ( kR ) -1 , in 
agreement with the interference concept of CB. Com- 
putations for significantly larger optical diameters will 
be needed to reach the (/< x transport mean free 
path) -1 scaling. 

2. Increasing absorptivity of the constituent particles can 
either enhance or suppress classical manifestations of 
the CB effect depending on the particle packing density 
and the real part of the refractive index. This duality 
illustrates again the extreme complexity of the effects 
of multiple scattering in many-particle groups (see, 
e.g., Ref. [46] where it is shown that absorbing Ray- 
leigh scatters can yield greater enhancement factors 
and deeper negative-polarization minima than non- 
absorbing Rayleigh scatters). 

3. In the case of very large values of the particle packing 
density, the scattering characteristics of the particulate 
volume start to exhibit behavior not predicted by the 
low-density theories of CB and RT. Nevertheless, the 
direct computer solutions of the Maxwell equations do 
demonstrate that the classical low-density predictions 
can survive (at least in a qualitative and even semi- 
quantitative sense) volume packing densities typical of 
particle suspensions and particulate surfaces encoun- 
tered in natural and man-made conditions. 

Furthermore, the accumulated body of evidence sup- 
ports the conclusion that CB is the likely cause of the 
uniquely narrow BOE and POE observed for a class of 
high-albedo Solar System objects [2-4]. Another impor- 
tant implication is that the azimuthal asymmetry of the 
polarized backscattering cone and the POE observed in 
the laboratory for densely packed particulate media 
[5,9,10] are indeed caused by the effect of CB. 


Finally, our numerically exact results pose a funda- 
mental question as to why the various manifestations of 
CB are so remarkably immune to packing density. The 
definite qualitative answer to this question may not be 
immediately obvious. However, one may speculate that 
even in densely packed discrete random media, the partial 
multiply scattered wavelets that involve widely separated 
particles still provide a significant cumulative contribu- 
tion to the total scattered signal and thereby make the 
classical multiple-scattering and CB features quite 
pronounced. 
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